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The dynamical properties at T = of the one- dimensional (ID) s — 1/2 nearest-neighbor (nn) 
XXZ model with an additional isotropic next-nearest-neighbor (nnn) coupling are investigated by 
means of the recursion method in combination with techniques of continued-fraction analysis. The 
focus is on the dynamic structure factors Szz{q,^^) and Sdd{<1,^), which describe (for q = n) the 
fluctuations of the Neel and dimer order parameters, respectively. We calculate (via weak-coupling 
continued-fraction analysis) the dependence on the exchange constants of the infrared exponent, the 
renormalized bandwidth of spinon excitations, and the spectral- weight distribution in SzziiTjU;) and 
SoDiir ,uj), all in the spin-fluid phase, which is realized for planar nn anisotropy and sufficiently 
weak nnn coupling. For some parameter values we find a discrete branch of excitations above the 
spinon continuum. They contribute to Szz{q,uj) but not to SDD{q,(^)- 
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I. INTRODUCTION 

Quantum many-body systems with competing interac- 
tions are apt to exhibit ordering tendencies in the ground 
state of which no trace exists in the presence of any 
one couphng alone and which are impossible to predict 
on the basis of a classical model with both couplings. 
Among the many different model systems where this 
phenomenon is manifest, the Heisenberg antiferromagnet 
with nearest-neighbor (nn) and next-nearest-neighbor 
(nnn) coupling on a bipartite lattice is a prominent ex- 
ample. For coupling strengths of comparable magnitude 
on the square lattice, various ordering tendencies includ- 
ing Neel order, coUinear order, dimer order, twist order, 
and chiral order are in competition with each other and 
with disordering tendencies .such as embodied by the res- 
onating valence-bond state.EJ 

The one-dimensional (ID) version of this model had 
gained prominence many years before, when the ground- 
state for one particular ratio of the coupJing constants 
was found to be a pure dimer state.BEl Spontaneous 
dimerization is a true quantum phenomenon. Sev- 
eral studies have succeeded in mapping out the zero- 
temperature phase diagram of ttia model in one or the 
other extended parameter spaceJjTJ albeit with little em- 
phasis on dynamical properties. In this paper we study 
the r = dynamics of the Hamiltonian 
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parametrized by the two coupling constants A = Jz/Ji_ 
and A = J2/ J±- Extending the parameter space to cases 
with uniaxial anisotropy in the nn coupling offers the 
advantage that the model with A — A_p 0, which reduces 
to a system of free lattice fermions,Elc can be used as a 
convenient starting point for weak-coupling approaches. 

A schematic representation of the T = phase diagram 
is shown in Fig. |^. The free-fermion point (A = 0, A = 0) 
is located in the middle of the spin-fluid phase. Here the 
system is a Luttinger liquid. The ground state is critical, 
i.e. the excitation spectrum is gapless, and the static 
spin and dimer correlation functions decay algebraically 
with exponents depending on the coupling constants of 
the two types of interaction. 

For A < the nnn coupling strengthens the cor- 
relations produced by the nn coupling and thus rein- 
forces the prevailing ordering tendency. Therefore, we 
expect to find ferromagnetic long-range order in the re- 
gion (A < 0, A < —1) and antiferromagnetic long-range 
order in the region (A < 0, A > 1). For — 1 < A < 1 the 
nn interaction alone is known not to support any kind 
of long-range order. No phase transition is suspected to 
occur if a negative nnn coupling is added. 

For A > the nn and nnn couplings are in competition 
with each other. The latter thus frustrates the ordering 



tendency of the former. One predictable consequence is 
that the boundary of the ferromagnetic phase, which is 
located at A = —1 for A < l/4_bends to the left in 
Fig. ^ for increasing nnn coupling.E5l Increasing amounts 
of uniaxial anisotropy are necessary to stabilize the spin 
alignment. A simple spin- wave stability criterion yields 
the expression 1/A = 4[|A| — \/A^ — 1] for the phase 
boundary in that region. 

In the region of planar anisotropy, |A| < 1, the in- 
terplay of the two competing forces in the presence of 
quantum fluctuations produces the dimer phase. Curi- 
ously, within the dimer phase therje jfixists one point in 
parameter space, (A = 1, A = l/2)p'0 where the ground 
state is a pure dimer state with no correlated fluctuations, 
notwithstanding the fact that this phase owes its very 
existence to quantum fluctuations. However, the T = 
dynamics at this point turns out to be only marginally 
simpler than elsewhere in the dimer phase, where j-the 
structure of the ground state is much more complexJlj 

In the region of uniaxial antiferromagnetic anisotropy, 
A > 1, there exists a Neel phase. Like the ferromagnetic 
phase, it is destabilized by sufRciently strong competing 
nnn coupling. The phase boundaries between the spin- 
fluid, dimer, and Neel ph|ases as sketched in Fig. |l| were 
flrst proposed by HaldaneEl based on a continuum fermion 
theory. The impact of the transitions between the Neel, 
dimer, and spin- fluid phases onlhe T = dynamics will 
be explored in a separate study.E3 



II. SPIN AND DIMER FLUCTUATIONS 

The Neel order parameter (staggered magnetizat ion ) 
and the dimer order parameter for the model system (1.1) 
are given, respectively, by the operators 
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Neel order in its purest form is realized in the states 

m) = \n'i ■■■ i) , 1*2) - uti • • • T) , (2.3) 

and pure dimer order in the states 



l*f> 



[l,2][3,4]...[iV- l,iV] 
[2,3][4,5]...[iV,l] . 



(2.4) 



where the singlets [I, I -|- 1] = {| ti) - | iT)}/\/2 are 
formed by pairs of nn spins. In the ground state of (pTl]), 
the former is realized at (A = 0, A = 00) and the latter 
at (A ^ 1/2, A = 1). 

Neel ordering manifests itself in the two-spin correla- 
tion function {SlS[_^_^) and dimer ordering in the four- 
spin correlation function {DiDi^n}- For translationally 
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invariant and orthonormal linear combinations of the 



s ym metry-breaking Neel states (2^) and dimer states 
(H, they are (5f 5f+„) = (AA+n) = (n ^ 

0), and (5f5f+„) = (|n| ^ 1), (AA+n)-(A)(A+„) = 
1)" {n ^ 0), respectively, and reflect the long-range 
nature of the two types of ordering. In the free-fermion 
ground-state at A = A = 0, by contrast, the same corre- 
lation functions decay algebraically,H 
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[cos(ri7r) — 1] 



(2.5) 



where Saa{q) = (^l^g^gl^") is the static structure factor 
(integrated intensity). 

For the reconstruction of S'aaC'Zjw) based on a num- 
ber of coefficients A^(g) extracted from the finite-size 
ground-state wave function |$), we employ techniques of 
continued-fraction analysis described previously in the 
context of other applications - one pertaining to the 
strong-coupling _efiginiec3 and the other to the weak- 
coupling regimeO of a given quantum many-body sys- 
tem. 



(AA+n) - (A)(A+„) 
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In this study we use the recursion methodt^HlJ to inves- 
tigate the dynamical (i.e. frequency-dependent) order- 
parameter fluctuations as probed by the dynamic spin 
structure factor S^zIIt^) a-nd the dynamic dimer struc- 
ture factor Sdd{<It^), i-e. by the function 



(2.7) 



where Aq stands for the spin and dimer fluctuation op- 
erators, 
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(2.8) 



For the calculation of the dynamic correlation function 
(<I>|Aj(t)Ag|$) in the ground state |$) of the system, it 
can be based on an orthogonal expansion of the-dynam- 
ical variable Aq{t) (Liouvillian representation)^^ or on 
an orthogonal expansion of tlm-.wave function Aq{~t)\^) 
(Hamiltonian representation) .E3 The algorithms of both 
representations, which have beeja, described and illus- 
trated in a recent monograph,EZI produce equivalent 
data. These data are expressible most concisely in terms 
of a sequence of continued- fraction coefficients A^(g), 
A^(g), . . . for the relaxation function 
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A^(g) 

... 



which is the Laplace transform of the symmetrized 
and normalized correlation function 5R(<I>|^J(t)Ag|<I')/ 
l<^\A\Aq\^). The r = dynamic structure factor (E 



is then obtained from (2.9) via 
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III. RESULTS AND INTERPRETATION 

For the special case A = A = of the spin model (|]|), 
dynamic correlation functions can be calculated exactly. 
The evaluation is particularly simple for the two dynamic 
structure factors of interest here, because the operators 
involved, Sq aiuLD^, are density operators in the fermion 



representation:E3E 



4J2 sin2((j/2)~w2 



4J2 sin2(g/2)- 
./isin2(q/2) 



(3.1) 



(3.2) 



for J_l| singi < w < 2J_l| sva{q/2)\. Both quantities have 
the same excitation spectrum — the particle-hole contin- 
uum of free fermions. For other dynamical variables such 
as Sq , the evaluation of dynamic correlation functions is 
much more involved, and the results have a much more 
complicated structure with dynamically relevant excita- 
tion spectra of unbounded support E3 

In the XXZ antiferromagnet with planar anisotropy 
(A = Q-O < A < 1), the spectral weight in Szziq,^^) is 
knowncJ to be dominated by a continuum of unbound 
two-spinon states with upper and lower boundaries 

eL(g) = JA|sinq| , ec;(q) = 2JA|sin(q/2)| , (3.3) 

where Ja = ttJ± smd/2'd, cosd — A. The A-dependence 
of Szz{Qt^) has already.,been investigated by several cal- 
culational techniques.Ea No equivalent results exist for 
Szzili^) iu the extended parameter space (A, A) or for 
Sdd{<1i^) anyw here in this par ameter space. Given the 
exact solutions (|3.l[) and ( |3.2| ), the case A = A = 
presents itself as a convenient starting point for a weak- 
coupling continued-fraction analysis of the-,two dynamic 
structure factors at |A| < 1 and |A| < 1.N 



A. Correlation exponent and spin velocity 

Haldane's continuum analysis! predicts that the expo- 
nents which characterize the power-law decay of the spin 



3 



and dimer correlation functions (5'fS'f^„) ^ (— l)"n , 
{DiDi+„) - {Di){Di+„) ~ (-l)"n-''^, are the same: 
Vz — Vd- a calculation taking into account both the 
backscattering and the umklapp terms in the interaction 
leads to a pair of scaling equations for the dependence of 
rjz on A and A. If the umklapp terms are neglected, the 
analysis yields the explicit result 



ry,(A, A) = 2 - 8A)/(^ + 4A) 



(3.4) 



This expression is expected to be most accurate for the 
special coupling ratio A/ A = 1/6, where the umklapp 
terms are absent in the continuum Hamiltonian. Indeed 



the high accuracy of (3.4) is demonstrated by the fact 
that its value rjz — 1.006 at the boundary to the Neel 
phase (A — 1,A = 1/6) misses the supposedly exact 
value = 1 by less than one percent. It is not clear how 
accurate expression ( |3.4[ ) is for coupling ratios A/A p4 
1/6. At A = it can be tested against the exact resulted 
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The two expressions agree to leading order in the cou- 
pling constant: -q^ = 2 ~ 4A/7r + ©(A^). 

An independent way of determining the correlation ex- 
ponent "qz for arbitrary coupling ratios A/ A in the region 
of very weak interaction is provided, as will be demon- 
strated in Sec. III.C, by the weak-coupling continued- 
fraction analysis of the infrared exponents in the dynamic 
spin and dimer structure factors. 
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The continuum analysis suggests that the two infrared 
exponents are identical and related to ry^: 



Pz^ Pd^ Vz 



(3.7) 



The renormalized spin velocity Vs{A,A) is another 
quantity for which the continuum analysis predicts an 
explicit result j3 



Vs{A,A) = J_i(H-2A/7r-4A/7r) 



(3.8) 



It can be checked against the renormalized bandwidth 
LUo{A.,A) of S'z2(7r,a;) or Sdd{''^,'^) as obtained (in 
Sec. III.C) via weak-coupling continued-fraction anal- 
ysis. Both quantities are expected to have the same de- 
pendence on the coupling constants except for constant 
factors. At A = we have luq — eu{TT) = 2Ja and 
Vs = [deuiq)/dq]q=o = ^A, hence 



UJn = 2w, 



(3.9) 



B. Weak-coupling continued-fraction coefficients 



The exact expressions (3.1) and ( p.2[ ) can be recovered 
for g = TT from the Afc-sequences 



(3.10) 



by direct evaluation of ( |2.9|) . These Ak values pertain to 
the infinite system. They are exactly reproduced up to 
k = N/2 - 1 for S'^z(7r,w) and up to fc = N/2 - 2 for 
Sdd{t^i cj) by the recursion method applied to a chain of 
iV sites with periodic boundary conditions. 

Weak couphng (|A| <C 1, |A| <C 1) produces system- 
atic d eviations of the A^'s from the reference sequences 
(pJQ). They arc illustrated in the four panels of Fig. 
p| for both dynamic structure factors and both types of 
interaction. In each case the A^'s for Ka < k < Kw ex- 
hibit a pseudo-asymptotic behavior. In panel (a) it starts 
at Ka — 2, in panel (b) at Ka = 3, and in panels (c), (d) 
at Ka = 1- The number Kw marks the beginning of the 
crossover from zero growth to power-law growth, which is 
most conspicuously observable in panel (d). Kw becomes 
smaller with increasing coupling strength. When we have 
Kw — Ka, we are in the strong-coupling regimcEj 

In the weak-coupling regime, the systematic deviations 
from (3.1C) are of two kinds: (i) The A^-sequence tends 
to converge toward a higher or lower value A^^ as k 
increases toward Kw- (ii) The A2k and the A2fe_i ap- 
proach the pseudo-asymptotic value A^'' from opposite 
sides. In Fig. || we observe that the direction of the 
shift, which changes with the sign of either interaction, 
is the same for both dynamic structure factors but oppo- 
site for the two types of coupling, whereas the alternating 
pattern, which also changes with the sign of either inter- 
action, has the same direction for both functions and for 
both types of coupling. 

For a quantitative analysis of the two pseudo- 
asymptotic effects we must select the string of A^ 's care- 
fully. In Sdd{t^,^), the value Kw is considerably smaller 
than in ^^^(Tr,^) for given coupling strength, which re- 
stricts the weak-coupling analysis of Sdd{t^, w) to a nar- 
rower range of couplings. 

The shift of the pseudo-asymptotic value A^'' [ef- 
fect (i)] describes the renormalized bandwidth loq of the 
dynamically predominan|t-2-particle continuum of lattice 
fermions via the relationllZI 



(3.11) 



For pure nn coupling (A — 0, |A| < 0.1), the band-edge 
frequency loq inferred from the average of A2, A5 was 
shown to reproduce the exact bandwidth eul^) = 2 Ja of 
the spinon continuum (3.3) very accuratelyO 

The alternating approach of the A^ 's to war d A^^-* [ef- 
fect (ii)] describes an infrared singularity ( |3.6| ) in the two 
dynamic structure factors. For a truly convergent Ak- 
scquence, the singularity exponent /3 is governed |by|-t|he 
leading term of the large-fc asymptotic expansion:ll3'ca 



Ak^Aoo[i-{-ir^ 



(3.12) 



In Ref. we proposed and tested two procedures (av- 
eraging and extrapolation) for extracting the exponent 
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from a finite Afc-sequence. For the nn case (A = 0, 
|A| < 0.05), a benchmark test showed that the expo- 
nent Pz inferred from the coefficients A2,...,A5 repro- 



duces [via (3.7)] the exact correlation function exponent 
(3.5) with reasonable accuracy. 

The patterns in Fig. ||make it clear that in the (A, A)- 
plane there are sectors with compressed (Acjq < 0) and 
expanded (Acjq > 0) bandwidth, and sectors with diver- 
gent (/? < 0) and cusp {(3 > 0) singularities. We have 
determined these boundaries between the resulting four 
sectors near the free-fermion point (A = 0, A = 0) by a 
systematic investigation of the A^-sequences along sev- 
eral lines in parameter space. 

One such series of sequences for 5*22(77, w) is displayed 
in Fig. 1^. Between A — —0.005 (bottom sequence) and 
A = 0.006 (top sequence) at fixed A — —0.001, we ob- 
serve a gradual reversal of both patterns (i) and (ii) . At 
A ~ —0.002, the average Afc goes from negative to posi- 
tive, implying a corresponding change in sign of Aujq . At 
A ~ 0.002, the average of A2jt_i — goes from positive 
to negative, implying a corresponding change in sign of 
the infrared exponent P^- The two sequences closest to 
the pattern changes are highlighted by full symbols. 

The lines in parameter space on which the interaction 
leaves the bandwidth Wg or the exponent (3z unchanged 
depend somewhat on the Afc-strings used for the analy- 
sis. This is illustrated in Fig. ^. The expectation is that 
the results improve as we shift the A^-string to higher ip^ 
dices fc, where non-asymptotic effects become weaker.E3 
Our best results, inferred from the string A4,...,A7, 
yield sector boundaries at coupling ratios A/A ~ 0.52 for 
Awq = and A/A ~ —0.44 for /3z = 0, in fair agreement 
with the corresponding sector boundaries A/ A = 0.5 for 
Avs = and A/A = —0.5 for Arjz — 0, res pect ively, 
predicted by the continuum results ( ^ ) and (3.4) and 
shown as long-dashed lines in Fig. 0. 



C. Infrared exponent and bandwidth 

Figure ^(a) shows the weak-coupling continued- 
fraction results for the infrared exponent f3z{A, A) at cou- 
pling strengths along the circle 



A = F cos 61 , A = F sin 6* 



(3.13) 



with radius F = 0.01 around the frcc-fcrmion point in 
the (ZY^)-plane. The three data sets from different A^- 
stringscZl fall onto sine-like curves in this representation, 
as do es the prediction inferred from the continuum result 
(3.4), shown here as solid line. 

The best overall agreement between our data and Hal- 
dane's result occurs at angles 9 ~ 9.5°, 189.5°, which 
correspond to the special coupling ratio A/ A = 1/6 with 
vanishing umklapp terms in the continuum Hamiltonian. 
However, the agreement between the dashed line and the 
data set from A4, . . . , A7, which is least affected by non- 
asymptotic coefficients, is remarkably good at all angles. 



The dependence of the renormalized bandwidth luq 
on the same angular variable 9 is shown in Fig. ^(b) 
in a comparative plot of three sets of weak-coupling 
cont inued- fraction data and the pred iction inferred via 
( |3.9D from the continuum result (3.S) for the renormal- 
ized spin velocity. Again the agreement improves as the 
Afc-string is shifted toward higher indices, where it is less 
affected by non-asymptotic coefficients. 

When the exponent and bandwidth are determined for 
a circle of much smaller radius, F = 0.0001, we find the 
same angular dependence of Pz and ujq and the same 
(rescaled) amplitude within a 5% margin of error. This 
confirms that the A^'s used in this analysis are free of 
weak/strong-coupling crossover effects. 

We have carried out the same analysis for the renor- 
malized bandwidth ujq and the infrared exponent Pd of 
Sooiq,^)- The results for parameter values on a circu- 
lar line with radius F — 0.0001 are displayed in Fig. |^. 
This radius had to be chosen much smaller than in the 
case of 5*22(77, cj) in order to ascertain that at least the 
coefficients A2, . . . , Ag are free of crossover effects. 

The bandwidth data from two different A^ strings de- 
picted in Fig. |^(b) are in near perfect agreement with 
each other and with the result inferred from (^|^) . Like- 
wise the data for the exponent Pd as shown in Fig. |^(a) 
exhibits the same characteristic oscillation as already ob- 
served in Fig. J^a) for Pz- The phase of the oscillation 
predicted by (3^) is accurately reproduced by our data. 
There are significant deviations in the amplitude, which 
are attributable to the strongly non-asymptotic nature of 
the coefficients A2 , A3 , but the trends indicated by the 
two data sets are in the right direction. 

The angular dependences of the exponent and band- 
width data shown in Figs. ^ and |^ are observed to be 
out of phase by some 55°, giving rise to the four sec- 
tors of weak-coupling dynamical behavior discussed in 

the context of Fig^ 0. This is consistent with Hal- 

dane's predictions ( ^.4| ) and ( |3.8D for the correlation ex- 
ponent and the renormalized Fermi velocity, respectively. 
These quantities have extreme values at tan 9 — 2 and 
tan^ = —2, respectively, which corresponds to angles 
63.4° and 116.6° = 63.4° -I- 53.2°. 

The characteristic oscillations of the infrared-exponent 
data are also expected to be present (with the same 
phase) in the angular dependence of the static spin and 
dimer structure factors at the critical wave number q = tt. 
In the free-f erm ion li mit, the Fourier transform of the ex- 
act results (2.5) and (2.6) yields: 



522(g) = SnD{q) = -[\q\-\sinq\]. (3.14) 

27r TT 

Both functions have a cusp-like maximum at g = 7r, 
which reflects the critical spin and dimer fluctuations, re- 
spectively, in the ground state. It is reasonable to expect 
that the variation of 522 (tt) and SooiT^) with the inter- 
action in the weak-coupling limit is synchronized with 
the variation of the correlation exponent ijz — rju, and 
hence with the infrared exponent Pz — Pd- 
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In Fig. we have plotted the 6'-dependence of S'zz(7r) 
and Sdd{t^) for fixed couphng strengths T = 0.0001 . 
Both quantities exhibit the characteristic sinusoidal be- 
havior. The amplitude is proportional to T in each quan- 
tity. Interestingly, the phases are different. Whereas 
SDoiT^) varies in phase with the infrared exponent, 
Szz{t^) does not. 



D. Reconstruction of ^^^(Tr,^;) and 51515(71,0;). 

The third major feature (in addition to the renormal- 
ized bandwidth and the infrared exponent) characterizing 
the line shape of the dynamic spin and dimer structure 
factors at the critical wave number g = tt is the detailed 
spectral-weight distribution near the band edge. For the 
pure nn case (A = 0), we found that A 7^ causes a re- 
distribution of the spectral weight in SzzIt^,^) near the 
band edge and (for A < only) thcj-appearance of a dis- 
crete spectral line outside the band.Ej This discrete state 
was identified (by Bethe ansatz) to belong to a branch of 
bound spin complexes. 



For the more general model (IT) with A 7^ 0, it is not 
known for which parameter values such discrete states 
exist. We have investigated this quest ion f or parame- 
ter values (A, A) on the circular line ( 3.13| ) by means 
of a weak-coupling continued-fraction reconstruction of 
Szzi"^, t^) based on the coefficients Ai, Ak^ and a ter- 
mination function which incorporates the renormalized 
bandwidth luq and the infrared exponent Pz- A detailed 
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description of the procedure can be found in Refs. ^ 

The dependence on the angular parameter 6 of the 
reconstructed 6*22 (7r,a;) is displayed in Fig. ^. As 9 
increases from zero to 60° [panel (a)], the bandwidth 
shrinks and the infrared divergence gains strength. In 
addition to these two features we observe the emergence 
at ~ 15° of a discrete spectral line above and out- 
side the band. In the reconstructed relaxation function 
the discrete state is represented by an isolated pole. For 
z = e — iLu with e = 0.0001 it has a nonzero width. At 
6 = 15°, the discrete state is barely distinguishable from 
the band edge. 

Between 9 = 60° and 9 = 120° the bandwidth con- 
tinues to shrink, while more and more spectral weight 
is transferred from the continuum to the discrete state. 
The infrared divergence, which has reached its maximum 
strength at ~ 60°, weakens over this parameter range. 
The further evolution of the line shapes near the band 
edge and at small frequencies is shown in panel (b). The 
infrared divergence disappears at ~ 150° and turns into 
a cusp singularity (see inset). The bandwidth of the con- 
tinuum expands while the discrete state moves closer to 
the band edge and slowly loses spectral weight. 

These trends continue between 6 = 195° and 6 = 240° 
as shown in panel (c) . The depletion of spectral weight 
at small uj becomes more and more prounouced, and the 
discrete state merges with the band edge at ~ 210°. 



Nothing dramatic happens to the line shape on the last 
stretch of the circle. The range of parameter values where 
a discrete state is observed, 30° < 6 < 210°, roughly co- 
incides with the range where the renormalized bandwidth 
is compressed. In the fermion representation, this is the 
region, where the Fermi velocity is renormalized down- 
ward. 

In the weak-coupling reconstruction of Sdd{'^,'^), 
which starts f rom t he expression (^^) and the constant 
Afc-se quence ( 3.10|) , the effects of a weak interaction 
( p.l3 ) on the line shape are similar in two aspects yet dif- 
ferent in a third aspect, (i) The low- frequency behavior 
is governed by a weak infrared singularity which switches 
from a divergence to a cusp in accordance with the expo- 
nent data presented in Fig. ||(a). (ii) The only noticeable 
change at the band edge is the variation of the contin- 
uum boundary in accordance with the bandwidth data 
presented in Fig. ^(b). (iii) However, no noticeable rear- 
rangement of spectral weight near the band edge occurs 
unlike what has been observed in 6*22 (7r,a;). In particu- 
lar, the discrete state found in SzziT^it^) does not carry 
any spectral weight in SuuinjUj). It does not play any 
role in the dimer fluctuations. 
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'^^ Obviously, the weak-coupling continued-fraction results for 
the infrared exponent Pz and the renormalized bandwidth 
ll>o from a string At^, A^^ of coefficients are vulnerable 
to systematic errors. The two main sources of systematic 
error are the following: (i) If ki is chosen too small {ki < 
Ka), then the first one or several Afc's are too strongly 
non- asymptotic in character and introduce a bias into the 
pseudo- asymptotic analysis, (ii) If fc/ is chosen too high 
(kf > Kw), then the last one or several A^'s are affected 
by the crossover to nonzero growth rate, which introduces 
a different kind of bias. 

In Ref. ^ we showed that the accuracy of the weak- 
coupling continued-fraction results for the case A = can 
yet be improved substantially if we use extrapolation in- 
stead of averaging. However, since the data do not lend 
themselves to extrapolation for all coupling ratios, we em- 
ploy the more robust averaging method for all cases here. 



FIG. 



1. Zero-temperature phase diagram of the model sys- 
tem (1.1) with coupling constants A — Jz/J± and A = J2I Ji_ 



in a schematic representation. 



FIG. 2. Continued-fraction coefBcients for the dynamic 
structure factors Szz{q,^jo) (left) and Sooiqjio) (right) at 
T = in the weak-coupling regime with A = (top) or A = 
(bottom). The energy unit is J_i = 1. 



FIG. 4. The spin-fluid phase in the weak-coupling regime 
is divided into four sectors with infrared exponents /3z and 
bandwidth renormalization Aujq of opposite sign. The solid 
lines denote the sector boundaries derived from the strings 
A2,...,A5 (circles), A3,...,A6 (squares), and A4,...,A7 (tri- 
angles). One solid line is extended short-dashed into the 
weak/strong-coupling crossover region of coefficient A7. The 
long-dashed lines with slope A/A — ±1/2 are the sector 
boundaries predicted by the continuum analysis. 



FIG. 5. (a) Infrared exponent /3z of (tt, tj) and (b) renor- 
malized bandwidth loq of the two-spinon continuum at differ- 
ent coupling ratios A/ A — tan 6 and fixed coupling strength 
VA2 + A2 = r = 0.01. The squares, circles, and triangles 
represent the results of a weak-coupling continued-fraction 
analysis based on three different Afe-strings. The solid lines 
in (a) and (b) are derived from the continuum results (3.4) 
and (3.8), respectively. The energy unit is Jx = 1. 



FIG. 6. (a) Infrared exponent /3d of SooiT^ji^) and (b) 
renormalized bandwidth of the two-spinon continuum at 
different coupling ratios A/ A = tan 9 and fixed coupling 
strength %/ A^ + A^ = F = 0.0001. The squares and circles 
represent the results of a weak-coupling continued-fraction 
analysis based on two different Afe-strings. The solid lines in 
(a) and (b) are derived from the continuum results (3.4) and 
(ETq), respectively. The energy unit is Jx = 1. 



FIG. 7. Integrated intensities 5'zz(7r) and ^^(Tr) of the 
spin and dimer order-parameter fluctuations, respectively at 
different coupling ratios A/ A = tan 8 and fixed coupling 
strength F = VA^ + A^ = 0.0001 as derived from the 
ground- state wave function for TV = 16. 



FIG. 8. Spin dynamic structure factor Szz (7r,tj) at dif- 
ferent coupling ratios A/ A — tan 9 and fixed coupling 
strength V A^ -I- A^ = F = 0.01, specifically for the angles 
(a) 9 = 0°,15°,...,60°, (b) 9 = 135°, 150°, . . . , 180°, (c) 
195°, 210°, . . . , 240°. The main plot depicts the line shapes 
near the band edge and the inset the line shapes at low 
frequencies. All curves have resulted from a weak-coupling 
continued-fraction reconstruction based on the coefficients 
Ai J. . . , A7 and a compact a-terminator as explained in Refs. 
Ilj^. The energy unit is Jx = 1. 



FIG. 3. Afc-sequences for the dynamic structure factor 
Szz{-K,uj) for A = -0.001, A = -0.005, -0.004, 0.006. 
The energy unit is Jx = 1. 
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